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ABSTRACT 

The Two degree Field Galaxy Redshift Survey (2dFGRS) is used in conjunction with 
the Two Micron All Sky Survey (2MASS) Extended Source Catalogue (XSC) to study 
the near-infrared light and stellar mass content of the local Universe. Mock galaxy 
catalogues, constructed from cosmological N-body simulations and semi-analytical 
galaxy formation models, are used to gauge the accuracy with which quantities 
can be recovered. The mean luminosity densities of the Universe are found to be 
pj = (3.57 ±0.11) X 10**/iLqMpc-3 and pKs = (7.04 ± 0.23) x lO^/iieMpc-^ (sta- 
tistical uncertainty only, and not accounting for the 2MASS low surface brightness 
incompleteness). Using the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) cat- 
alogue, the group mass-to-light ratio in the Ks band is found to increase by a factor 
of ~ 3 when going from groups with total 6j-band luminosities of 3 x 10^" 
to the richest clusters. These clusters have typical dynamical mass-to-light ratios of 
T K ~ SO/iTq. Galaxy luminosities are used to estimate stellar masses. Taking into ac- 
count the bias introduced by uncertainties in estimating galaxy stellar masses, a value 
of fJstars^ = (0.99 ±0.03) X 10"^ is measured, assuming that a Kennicutt stellar initial 
mass function (IMF) is applicable to all galaxies. Changing this to a Salpeter stellar 
IMF gives ristars/i ~ 2.1 X 10"'^. The 2PIGGs are then used to study the distribution 
of the stellar content of the local Universe among groups of different size. The three 
main conclusions arc: 

(1) a slowly rising stellar mass-to- i^s band light ratio is found with the clusters having 
the largest value of ~ 0.6Tq, 

(2) in contrast, the fraction of mass in stars decreases with increasing group size, 
reaching ~ 5 x 10~^/i for the rich clusters, and 

(3) in answer to the question posed in the title, most stellar mass is contained in Local 
Group-sized objects (M ^ 2 x 10^^ /i^^ M©) with only ^ 2 per cent in clusters with 

5 X 10"/i-1Mq. 

Key words: galaxies: groups - galaxies: haloes - galaxies: clusters: general - large- 
scale structure of Universe. 



1 INTRODUCTION 

Understanding the impact that environment has on a 
galaxy's ability to form stars is one of the main goals of 
cosmology. A very pertinent observational constraint is the 
manner in which stellar mass is distributed between galaxy 
systems of different size. This paper describes a measure- 
ment of this constraint, using a combination of the Two 
Micron AH Sky Survey (2MASS) Extended Source Cata- 
logue (Jarrett et al. 2000) and the 2dFGRS (CoUess et 
al. 2001, 2003) Percolation-Inferred Galaxy Group (2PIGG) 
catalogue (Eke et al. 2004). 



By-products of this quest, which are also reported here, 
include a recalculation of the local galaxy near-IR luminosity 
and stellar mass functions with the largest sample of galaxies 
yet used for this purpose. Mock galaxy catalogues are used 
in order to help quantify systematic uncertainties and biases. 
In addition, the dependence on group size of the stellar mass 
fraction and of the group mass-to-light and stellar mass-to- 
light ratios are studied and compared with previous results. 

Section 2 contains a description of the data used in this 
analysis. The near-IR luminosity functions and mean lumi- 
nosity densities are calculated in Section 3 and compared 
with recent determinations by Cole et al. (2001; hereafter 
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COl), Kochanek et al. (2001) and Bell et al. (2003). Group 
total mass-to-light ratios in the Ks band are reported in 
Section 4. These results extend the range of group sizes over 
which near-IR meiss-to-light ratios have been measured, rel- 
ative to the studies of Kochanek ct al. (2003), Lin, Mohr & 
Stanford (2003; hereafter LMS03) and Rines et al. (2004). 
Section 5 contains a description of the procedure used to 
calculate stellar masses given a set of galaxy luminosities in 
different wavebands. The accuracy of these measurements is 
gauged using mock galaxy catalogues produced from a com- 
bination of a dark matter N-body simulation and a semi- 
analytical galaxy formation model. This information is then 
used in Section 6, where the total galaxy stellar mass func- 
tion and fistars are estimated. 

The stellar content of the galaxy population is split by 
group size in Section 7, where group stellar mass-to-light 
ratios and stellar mass fractions are presented and com- 
pared with those found by LMS03. Finally, Section 8 pro- 
vides the constraint alluded to at the start of the Introduc- 
tion. Throughout this paper, it is assumed that Qm = 0.3, 

= 0.7 and the Hubble constant is written as Ho = 
100/ikms-^Mpc-^ 



2 OBSERVATIONAL DATA 
2.1 2dFGRS 

The two contiguous patches of sky covered by the 2dFGRS 
contain approximately 190 000 galaxies with a median red- 
shift of 0.11. For the purposes of the work performed in this 
paper, only the ~ 109 000 galaxies a.t z < 0.12 were used. 
At higher redshifts, less than half of the 6j-band luminos- 
ity of typical groups is contained within galaxies above the 
survey flux limit. The parent catalogue for the 2dFGRS has 
a very well quantified completeness. Norberg et al. (2002) 
and Cross et al. (2004) estimate that only ~ 9 per cent of 
the 6j-band luminosity is missing, while Driver et al. (2004) 
conclude from the MGC that the 2dFGRS-infcrrcd luminos- 
ity density is 3 per cent low due to missing low surface 
brightness galaxies (at all luminosities). 

In Sections 4, 7 and 8, where the galaxies are split by 
host group size, the 2PIGG catalogue (Eke et al. 2004) has 
been used to define the galaxy groups and their properties. 
Note that in Section 8, where the distribution of stellar mass 
throughout the group population is studied, single galaxies 
are also defined as groups. While such groups do not have 
a dynamically estimated mass, they still have a well defined 
total group luminosity. 



2.2 2MASS 

The near-IR data wore obtained from the 2MASS All-Sky 
Extended Source Catalogue (Jarrett et al. 2000), using 
the database available at http://irsa.ipac.caltech.edu/cgi- 
bin/Gator/nph-dd. A total of 43 553 z < 0.28 2MASS galax- 
ies were matched with those in the two contiguous patches 
of the 2dFGRS. The median redshift of the 2dFGRS galax- 
ies with 2MASS matches is ~ 0.095, slightly lower than the 
median of all 2dFGRS galaxies (~ 0.11). 

The total galaxy J-band magnitude was estimated from 



slope l,003±0,0aG 

14 - 




K (Loveday Kron) 

Figure 1. Comparison of the 2MASS A'-band Kron magnitudes 
with those measured by Loveday (2000). A straight hnc traces 
the least-squares fit to the data. The inset histogram shows the 
distribution of magnitude differences between the 2MASS and 
Loveday measurements, which is interpreted as the distribution 
of errors in the 2MASS values. The best-fitting Gaussian is also 
shown, and its parameters are given in the legend. 

the J-band Kron measurement (corrected for Galactic ex- 
tinction), including an offset of — O.lmag to account for flux 
outside the Kron aperture. This is approximate because 
of the variety of different surface brightness distributions 
among the galaxies, but it is a reasonable typical correc- 
tion (COl). The total galaxy Ks-ha.nd value is then inferred 
using the total J-band magnitude and the J — Ks colour 
within the 20 mag per square arcsec isophotal fiducial el- 
liptical aperture. Following COl, the -ffg-band Kron mag- 
nitudes (without the —O.lmag shift) were compared with 
the deeper pointed observations of Loveday (2000). Fig. 1 
shows that the 2MASS Ks-hand Kron values inferred from 
the J-band Kron magnitude and the isophotal colour are 
essentially unbiased, and have a scatter about the Loveday 
measurements that depends strongly on magnitude. This 
represents an improvement over the first 2MASS data re- 
lease used by COl, who needed to apply an offset to match 
the Loveday data. For the whole sample of 80 galaxies, the 
distribution of errors, shown in the inset histogram, has a 
best-fitting Gaussian standard deviation of 0.125mag. 

When creating mock catalogues with near-IR galaxy 
magnitudes, it is necessary to include the observational er- 
rors in order to determine the systematic bias that they 
introduce. The medians of the binned errors on the J-band 
total magnitude and the J — Ks colour are shown in Fig. 2 as 
functions of the total J-band magnitude. These have been 
constructed using the 1 a uncertainties returned from the 
2MASS database. As was apparent in Fig. 1, the observa- 
tional errors are greater for the fainter galaxies. 

The k + e corrections adopted for the J and Ks bands 

are 
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Figure 2. The errors on the 2MASS photometry are shown as 
functions of galaxy J-band total magnitude. Crosses and filled 
squares represent the medians of the J-band and Ks — J Icr 
uncertainties in each bin, i.e. the typical errors. The individual 
galaxy errors used to construct these medians are those contained 
in the 2MASS database. The two lines trace the quadratic fits 
given in the panel. 



k + e. 



and 



k + e. 



(g + 4g^) 
' 1 + 9^;2-8 



' 1 + 922-5 



(2.1) 



(2.2) 



respectively. These are fits to the mean of the k + e correc- 
tions computed for individual galaxies using stellar popula- 
tion models. A detailed description of how the k + e correc- 
tions are calculated can be found in COl. In brief, look-up 
tables oi bj, rp, J and K magnitudes as a function of red- 
shift were generated for a grid of star formation histories, 
using the stellar population spectral energy distributions 
provided by Bruzual & Chariot (2003). Observed galaxies 
were matched with theoretical star formation histories, thus 
allowing a k + e correction to be assigned to each galaxy 
using the theoretical model. These simple models are used 
in Section 5 to assign stellar masses to galaxies. 

Note that the mean V/Vmax = 0.513 for the 16 922 
galaxies with 11.8 < J < 14.6 and z < 0.12. In this ex- 
pression, V represents the volume surveyed out to the red- 
shift of each particular galaxy, and Vmax is the maximum 
volume in which each galaxy could have been detected. The 
isTs-band sample of 15 664 galaxies with 10.8 < Ks < 13.4 
and z < 0.12 has a mean V{zi)/V{znia:>i,i) = 0.517. These 
values suggest that the galaxy samples are not significantly 
incomplete (which would lead to values lower than 0.5) when 
these flux limits arc applied. The top panels in Fig. 3 show 
the full distributions of V{zi)/V{zinax,i) for the samples used 
in the two bands, and the lower panels show the differential 
number counts with the Euclidean slope removed for all of 
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Figure 3. The top row shows the distributions of V/Vmax for 
the J- (left) and Ks-ha.nd (right panel) samples. The lower row 
contains plots of the total number counts (AT) as a function of 
magnitude for these same two bands. Vertical lines illustrate the 
flux limits applied to the survey in this paper. 



the matched galaxies. These suggest that there is no strong 
systematic depletion of the galaxies near the adopted flux 
limits. Andreon (2002) and Bell et al. (2003) find a small 
incompleteness in the 2MASS survey arising from a deficit 
of faint, low surface brightness galaxies. While this makes 
a significant difference to the luminosity and stellar mass 
functions for the smaller galEixies, it has only a ~ 20 per 
cent impact upon the total luminosity or stellar mass den- 
sities inferred from the data. This will be discussed further 
in Section 6. 



3 LUMINOSITY FUNCTIONS 

With the galaxy total J- and Ks-hand magnitudes defined 
as described in the previous section, it is now possible to cal- 
culate the galaxy luminosity function in these bands. This 
analysis extends that by Kochanek et al. (2001) and COl 
by using over twice as many galaxies, and also by adopt- 
ing a different treatment of the statistical errors. In addi- 
tion, mock 2dFGRS catalogues have been analysed to help 
determine the accuracy to which the luminosity functions 
can be recovered. The procedure for creating these mock 
catalogues was described in Eke et al. (2004a), and basi- 
cally relies on populating dark matter haloes formed in the 
GIF N-body simulation (Jenkins et al. 1998) using the semi- 
analytical galaxy formation model described by Cole et al. 
(2000). In order to match the 2dFGRS 6j band luminosity 
function a luminosity-dependent colour-preserving shift is 
applied to each galaxy. However, the resulting model does 
not match the near-IR galaxy luminosity function. Thus, for 
the purposes of this paper, where the mock catalogue is be- 
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Table 1. Best-fitting Schechter function parameters for the different galaxy luminosity functions. The normalisations come from a 
fit to the 1/Vinax estimate of the abundances of the most luminous galaxies. Also listed are the adopted solar absolute magnitudes in 
the two bands and the mean Universal luminosity densities implied by the SWML functions. Only statistical uncertainties are quoted. 
The larger systematic uncertainties are discussed in the text. No account has been taken of the systematic impaxit of the errors in the 
2MASS magnitudes. However, this effect is smaller than the quoted statistical uncertainties. 
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Figure 4. J-band luminosity functions for mock catalogues and 
2dFGRS/2MASS data. The solid histogram traces the model 
galaxy J-band luminosity function, whereas the dashed line shows 
the luminosity function actually recovered from the mock cata- 
logue using a SWML estimator, normalised to agree with the 
l/Vmax luminosity function for 10.1 < logiQ[Lj/{h~^ Lq)] < 
11.2. Results for the SWML estimate from the real data are shown 
with the points, and the error bars arc inferred using a jack-knife 
method. The dotted line traces the best-fitting Schechter function 
derived from a fit to the SWML estimate over the same range 
of luminosities used above. The Schechter function parameters 
are given in Table 1. 



ing used to quantify errors in the recovered quantities rather 
than testing the galaxy formation model itself, an additional 
colour-preserving shift is applied when considering near-IR 
luminosities or stellar masses. In effect, each galaxy is made 
15 per cent brighter, this being the factor necessary to in- 
crease the mean Ks band luminosity density in the model to 
agree with that in the real data. This simulation cube filled 
with galaxies, hereafter referred to as the model, is then 
replicated and observed with suitable flux limits, masks etc. 
(Norberg et al. 2002), creating a mock catalogue. While the 
mean Ks band luminosity density is the same in mock and 
real catalogues, the luminosity function is not, so the near- 
IR flux limits applied to the mock catalogue are chosen to 



Figure 5. The equivalent of Fig. 4 for the Xg-band lumi- 
nosity functions. The SWML estimator has been normalised 
by matching to the l/Vniax luminosity function at 10.4 < 
logioliifs/C^-'Lo)] < 11.4. 



recover the same number of galaxies as are being considered 
from the 2MASS sample. 

Figs 4 and 5 show the J- and Ks-haxid luminos- 
ity functions derived from both the mock catalogue and 
the 2dFGRS-2MASS data. Only the results obtained us- 
ing a stepwise maximum likelihood estimator (SWML; Efs- 
tathiou, Ellis & Peterson 1988) are shown. A l/Vmax estima- 
tor yields lower numbers of less luminous galaxies because 
both mock and real observers happen to live in underdense 
regions of universe. However, this estimator does give a re- 
liable abundance at higher luminosities, where it probes a 
larger volume, and has thus been used to fix the normalisa- 
tion of the SWML curves. This was done using a minimi- 
sation in either of the intervals 10.1 < logiglLj /{h~^ Lq)] < 
11.2 or 10.4 < logio[iKs/(/i~^i0)] < 11-4. Errors on the 
recovered SWML luminosity functions arc dotorminod from 
a jack-knife procedure. This entails splitting the galaxy sam- 
ple, by right ascension, into ten subsamplcs containing equal 
surveyed areas, recalculating ten different estimates (ni) of 
the luminosity function by rejecting each subsample in turn, 
then defining the uncertainty as 
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Table 2. The near-IR luminosity functions recovered with the 
real data. Luminosities axe measured in h^^ Lq and the quoted 
abundances are dn/dlogjQL in (h~^Mpc)~^ . 
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where (n) represents the luminosity function obtained using 
the full sample. The J and Ks luminosity functions esti- 
mated from the mock catalogue provide an accurate recov- 
ery of the underlying model function (shown with solid his- 
tograms), although they do differ somewhat from the near- 
IR luminosity functions estimated for the real data, despite 
having the same mean luminosity densities. The deficit of 
high- luminosity galaxies in the mock catalogue represents an 
interesting discrepancy between the semi-analytical model 
assumed in the mock catalogue and the real Universe. 

Having normalised the SWML luminosity function us- 
ing the luminous end of the 1/Vmax estimate, a fit can 
be performed to determine the best-fitting Schechter (1976) 
function: 

0(L)dL = 0, ( A) %xp (- £) ^ . (3.2) 

In principle, the range over which this fit is performed need 
not be restricted to the luminous galaxies. However, the im- 
pact of any potential low surface brightness incompleteness 
in 2MASS (Andreon 2002) is likely to be apparent at lower 
luminosities, so the fit is restricted to the same range of 
luminosities used to normalise the SWML function. Thus, 
the derived parameters, which are listed in Table 1, are also 
those that would have been found by fitting to the 1/V^ax- 
estimate of the luminosity function. 



The mean Universal luminosity densities implied by the 
SWML-estimated luminosity functions are also included in 
Table 1. Statistical uncertainties come from the dispersion 
in the fits to the jack-knife samples. Note that, compared 
with the analysis of COl, the uncertainties on the luminosity 
function shape parameters (L* , a) are larger, because of the 
limited range of luminosities over which the fit is performed. 
In contrast, the statistical uncertainty on the normalisation 
is lower because of the increased sample size (both larger 
area and a deeper survey limit) and the choice of normali- 
sation method. Matching the abundance of high luminosity 
galaxies is a better way of reducing the variance than fitting 
to the counts, as COl did. 

The statistical uncertainties arc small compared to the 
systematic uncertainties. These include imperfections in the 
sample, due to incompleteness and misclassification in the 
parent 2dFGRS (Cross et al. 2004) and missing low surface 
brightness galaxies in the 2MASS (Andreon 2002; Bell et 
al. 2003). It is this last effect that is most important. Both 
Andreon (2002) and Bell et al. estimate that the 2MASS 
XSC misses ~ 25 per cent of the mean luminosity density 
in the Ks band. The mean near-IR luminosity densities in 
Table 1 are consistent with those of COl and Bell et al. 
(2003), but have smaller statistical uncertainty. 



4 ifs-BAND GROUP MASS-TO-LIGHT 
RATIOS 

Following Eke et al. (2004b), who studied the bj- and rp- 
band galactic content of groups of different size, the near-IR 
group mass-to-light ratio {Tk) is now considered. As the 
total galaxy stellar mass correlates more strongly with the 
near-IR flux than it does with the optical flux, Tk more 
effectively reflects the efiiciency with which stars have been 
formed in regions that end up in these different sized haloes. 



4.1 Results for the model 

Fig. 6 shows how the star formation efficiency varies with 
halo size. The model used to populate the mock catalogues 
shows a clear minimum in mass-to-light ratio at group sizes 
corresponding to about one L* galaxy, as was seen in the op- 
tical bands (Eke et al. 2004b), as well as a plateau for clus- 
ters with logjQ(Li,j/ /i^^ L0) > 11.5. The process of iden- 
tifying groups and inferring their mass-to-light ratios using 
a dynamical mass estimator inevitably introduces errors. A 
dotted line traces the mass-to-light ratio variation recovered 
from the mock catalogue. In selecting the groups to be used 
to make this curve, only local and isolated ones are taken, to 
avoid some sources of contamination. As a consequence of 
these restrictions, which are detailed in the caption of Fig. 6, 
the group mass-to-light ratio can be recovered to within 
~ 50 per cent for groups with log^Q(I/6j/ Lq) > 10.5. 
The main source of bias is due to interloping galaxies in the 
groups, which tend to increase the luminosity more than 
the mass (the groups are found in redshift space so inter- 
lopers have velocities similar to those of true group mem- 
bers). An increase of Tk by a factor of 3 in going from 
logjQ(L6j/ Lq) = 10.5 to 11.3 is recovered, which is com- 
parable to that present in the model haloes, albeit shifted 
to slightly larger halo luminosity. 
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Figure 6. The variation of the mock and 2PIGG Xg-band mass- 
to-light ratios with halo size, parametrized by the total fej-band 
luminosity. A horizontal line shows the mean mass-to-light ratio of 
the model, while the solid curve traces the variation of the binned 
median mass-to-light ratio in the model. The accuracy with which 
this is recovered is shown by the dotted line, and the points 
display the results obtained from the real 2PIGG sample. Only 
groups with no neighbours within a distance of dmin/ Mpc = 
2 -f [ 10 — logj^()(Li,j / hr"^ Lq)] and out to a maximum redshift of 
2,-,jjjx = 0.04 + 0.03 [logj^„(Li,j/ h^'^ Lq) - 10] arc considered. Er- 
ror bars show the 16th and 84th percentiles in each bin. A dashed 
line traces the results of Lin et al. (2003), mapped onto these axes 
using the median mass-to-6j light found by Eke et al. (2004b), and 
assuming that the clusters have a concentration of 5. Taking the 
median group total mass-to-light ratio of Eke et al. (2004b) maps 
the following luminosities logiQ[L(,j/(/i~^L0)] = (10,11,12) to 
masses logio[Mo/(h-iM"©)] = (12.0, 13.6, 14.7). 



4.2 Observational data 

The points in Fig. 6 represent the result of applying the 
same groupfinding, mass-to-light measuring, and group se- 
lecting procedure to the 2dFGRS as was applied to the mock 
catalogue. Again, a factor of ~ 3 increase in T is visible 
when going from the smaller groups to rich clusters. All 
of this trend occurs at logiQ(L6j/ Lq) < 11.3. Above 
logjQ(Li,j/ Lq) « 11.3, the ncar-IR group mass-to-light 
ratio remains approximately constant. This higher luminos- 
ity range corresponds to the range of cluster sizes in the sam- 
ple of Lin et al. (2003, LMS03). They find that the Ks band 
T500 oc MgQQ^*""®. This slope is steeper than suggested by 
the data presented here. These are in better agreement with 
the results of Kochanek et al. (2003) who find this expo- 
nent to be 0.10±0.09, despite going down to smaller groups 
(logiQ(Lj,j/ Lq) « 10.8) where the trend starts to be- 
come apparent in the 2PIGG results. 

For the 110 groups with log]^g(L6j / /i"'^ Lq) > 11.3, 
the median Tk = 89 ± 3 /iTq (statistical uncertainty 
only). If the 2PIGG mass-to-light ratios for such clusters 
are overestimated by the same factor as in the mock cat- 



alogues, then one should correct the median mass-to-light 
ratio downwards accordingly to yield Tx = 77 ± 3 /iTq. 
This is very similar to the value found by LMS03. For a 
sample of 19 clusters with kTx > 3.7 keV, they obtained 
Tx = 76 ± 4 /iTq using the 2MASS ii's-band group lumi- 
nosities with a statistical background correction, and masses 
inferred from X-ray data and the assumption of hydrostatic 
equilibrium. The LMS03 clusters extend to blue luminosities 
of logj„(L6,/ /i^^ Lq) ~ 13, off the right hand side of Fig. 6, 
and it is these higher luminosity clusters that bring the me- 
dian mass-to-light ratio up into agreement with the 2PIGG 
results. Kochanek et al. (2003) find Tx = 116 ± 46 /iTq, 
which is consistent within the statistical uncertainty. 

There are a number of potential systematic differences 
that blur the comparison between the cluster Jfg-band mass- 
to-light ratios reported in different papers, at a level of ~ 10 
per cent. These include the possibility of galaxies not trac- 
ing dark matter in real haloes as they do in the model 
universe and incompleteness in the parent 2dFGRS. The 
LMS03 value is measured at a radius, rsoo, enclosing a mean 
density of 500 times the critical value, whereas the group 
definition used here is tuned to the halo found by apply- 
ing a friends-of-friends group-finder, with a linking length 
equal to 0.2 times the mean interparticle separation, to a 
dark matter simulation. This algorithm finds groups out 
to a radius rvir t'2oo > fsoo- Rines et al. (2004) find 
that the variation of mass-to-light ratio with radius is slight 
in their sample of 9 rich clusters. They use mass profiles 
determined from the infall pattern of galcixies. Combining 
this with an assumption that the dark matter haloes have 
Navarro, Frenk & White (NFW, 1997) profiles with a con- 
centration of 5, they find that the mean mass-to-light ratios 
are T(< rsoo) = 78 ± 7 /iTq and T(< 7-200) = 88 ± 9 /iTq. 

4.3 Estimating fim 

Given the ATg-band mass-to-light ratio of clusters, an as- 
sumption as to how this compares to the mean Universal 
mass-to-light ratio and a measurement of the mean lumi- 
nosity density, yields an estimate of f2m • If real clusters, like 
the model clusters, have mass-to-light ratios that underes- 
timate the mean Universal value, then this exercise gives 
Qm = 0.21 ± 0.01 (statistical). This value is obtained after 
applying two largely offsetting ~ 10 per cent corrections, 
one assuming that the cluster mass-to-light ratio is overesti- 
mated and the other that it provides an underestimate of the 
universal value, so the statistical error is likely to bear little 
relation to the true uncertainty on this estimate. It is inter- 
esting to note that this estimate of S7ni is somewhat lower 
than that produced by Eke et al. (2004b) using a similar 
method in the 6j and rp bands (f2m ~ 0.27). 

While the total Ks band luminosity density in the mock 
is the same eis that in the real world, the mock clusters have 
mass-to-light ratios that are ~ 30 per cent greater than those 
in the real world. If the mock catalogues had more faithfully 
represented the real shape of the T^s-band galaxy popu- 
lation, then not only would more luminous galaxies have 
existed in the model but, perhaps, they would have prefer- 
entially occurred in rich clusters. The net effect would have 
been to make the cluster mass-to-light ratio a smaller frac- 
tion of the universal mean. Consequently a larger correction 
upwards would have been required, and when applied to the 
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real data a higher fim would have resulted. Given the sys- 
tematic difTcrences between the mock catalogue and the real 
data in the Ks band, it would be premature to say that this 
low apparent value of fim represents anything more than a 
challenge for semi-analytical models to place more i^Ts-band 
luminosity into clusters. 



5 STELLAR MASS ESTIMATES 

The stellar mass of a galaxy is estimated from its broad band 
magnitudes using simple stellar population synthesis mod- 
els, as described at length by COl. In summary, a grid of 
star formation histories is constructed, in which both the 
e-folding time, r, which governs the star formation rate, 
{tp oc exp(— f/r)), and the metallicity of the stars are varied. 
When combined with the assumption of a Kcnnicutt (1983) 
stellar initial mass function (IMF) and the stellar population 
synthesis models of Bruzual & Chariot (2003), look-up ta- 
bles of broad band magnitudes versus rcdshift can be gener- 
ated. Each observed set of galaxy magnitudes is matched up 
with the theoretical model that best reproduces these mag- 
nitudes at the redshift of the galaxy. The theoretical track 
can then be used to estimate the k+e correction that should 
be applied to calculate the absolute magnitude of the galaxy 
(see Section 2.2) and also to infer its mass-to-light ratio at 
z = 0. There are two major differences between the stellar 
masses computed here and those inferred by COl. First rr 
magnitudes have now been measured for the 2dFGRS ob- 
jects. Second, models are selected from the theoretical grid 
to produce unique predictions for colour versus redshift. This 
means, in practice, a reduced dynamic range of the grid in 
e-folding time and metallicity compared with that used by 
COl. The stellar metallicities used here are now closer to 
solar. 

With the availability of rp magnitudes, there is now 
more information to help assess which theoretical model best 
matches the colours of each observed galax:y. Throughout the 
rest of this paper, two different stellar mass estimators will 
be referred to. The J-band stellar mass will be that coming 
from using all four fluxes (6j, rp, J and Ks) to choose a 
model star formation history. The second estimate, the vf- 
band inferred stellar mass, will be identical to the J-band 
estimate when near-IR fluxes are available, but will also in- 
clude galaxies too faint to make the near-IR flux limit of the 
sample. In these cases the rp-inferred mass will be based on 
just the bj and rp data. These two stellar mass estimates 
should thus be considered to represent two different subsets 
of galaxies, with the rp-inferred stellar masses providing a 
deeper sample that reaches higher redshifts and intrinsically 
less massive galaxies than the J-band sample. 

Note that the choice of stellar IMF introduces a large 
systematic uncertainty in the recovered stellar masses. For 
instance, as COl showed, a Salpeter IMF assigns almost 
twice as much stellar mass per unit light as is inferred using 
the Kennicutt IMF that is used throughout this paper. It 
should also be stressed that the stellar masses referred to in 
this paper are those currently locked up in stars. This differs 
from the total mass of star formation by a factor that de- 
pends on how much stellar mass is recycled. Following Cole 
et al. (2000), a recycled fraction oi R = 0.42 is adopted 
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Figure 7. Errors in the inferred stellar masses using the rp (up- 
per panel) and ,/ band (lower panel) estimators. Tlie points show 
the median bimied accuracy with which stellar masses are recov- 
ered from the mock catalogue using the rp-band stellar masses. 
Error bars display the 16th and 84th percentiles. All galaxies at 
z < 0.12 have been used to calculate these quantities. The solid 
lines represent the mean ±lcr errors used to model this bias in 
Section 6. Dotted vertical lines delimit the region containing the 
central 80 per cent of the stellaj: mass (i.e. only 10 per cent comes 
from both regions not between dotted lines). 



throughout this paper to approximate this effect for a Ken- 
nicutt IMF. 

Another uncertainty in the estimated stellar masses is 
the appropriateness of the assumed exponential star forma- 
tion history. It is possible, however, to estimate the size of 
this error using the mock catalogues. Fig. 7 shows how well 
the recovery works for the rp- and J-band stellar masses 
referred to above, as a function of the true stellar mass in 
the model galaxy. The inclusion of magnitude errors con- 
tributes little to the range of stellar mass errors, which are 
instead dominated by the choice of the grid of star formation 
histories. While the grid of models assumes that the stars 
formed with a single exponential timescale, the more real- 
istic semi-analytical galaxies undergo much more complex, 
bursty histories. Consequently, some non-negligible error in 
the estimated stellar masses is to be expected. Rough fits to 
the typical errors in the stellar mass recovery and the scat- 
ter about this error are shown with the lines in Fig. 7. The 
equations of these lines are, for the rp band: 

J/ = 1.34 - 0.12a; (5.1) 

and 
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ay = 0.4 - 0.025a;, (5.2) 

where y = logig[Mstars(estimated)/Matars(true)] and x = 
logiQ[MatarB(true)/( M0)], and for the J band 

y = 1.18 - 0.10a; (5.3) 

and the same cr„ as above. 



6 THE TOTAL STELLAR MASS FUNCTION 

For each of the stellar mass estimates, it is possible to deter- 
mine the total stellar mass function in galaxies. The mock 
catalogues provide guidance on the size of possible system- 
atic errors arising from the measurement procedure. Both 
1/Vmax and bivariate stepwise maximum likelihood (SWML; 
Sodre & Lahav 1993; Loveday 2000) methods have been em- 
ployed. The former method is sensitive to large scale fluctu- 
ations in the galaxy density field and is thus less reliable for 
the low stellar masses that arc only seen locally. However, 
the high mass end with 10.4 < logio(Mstara/ h~'^ M©) < 11.4 
provides a good normalisation of the SWML method, which 
is insensitive to these density fluctuations but only returns 
a shape for the stellar mass function. A bivariate SWML 
method is necessary, with the variables being the inferred 
stellar mass and the gala^xy luminosity in the waveband 
defining the survey (i.e. 6j or J for the tf- and J-band stellar 
masses respectively). 

6.1 Systematic errors in the recovery 

Fig. 8 shows how well the stellar mass function can be recov- 
ered from the mock catalogue. Only galaxies with z < 0.12 
are used in constructing these functions. It is immediately 
apparent that the effect of the measurement errors in the 
galaxy stellar masses is to produce an overestimate of the 
abundance of galaxies at the higher stellar masses. This, in 
turn, is reflected in an overestimation of the contribution 
to the cosmic density from stars residing in galaxies. While 
the model value is Qstmsh = 1.12 x 10~^ (from integrat- 
ing under the solid black line multiplied by Mstars), and the 
fraction of universe that is actually probed by the mock sur- 
vey only has f2atara/i = 1-03 x 10~^ (dotted line) because of 
variations due to large scale structure, the recovered SWML 
values are Ostara/i = 1.46 x 10~^ (short-dashed line) and 
f^atara/i = 1.45 X 10~^ (loug-dashcd line) for the rp- and 
J-band measurements respectively. To give an idea of the 
size of the error that can be introduced by using the 1 /Vmax 
stellar mass function, rather than the suitably normalised 
SWML one, the J-band 1/Knax stellar mass function yields 
a value of Qstaieh = 1.26 x 10~^, ~ 15 per cent down on 
the SWML estimate, which is unaifected by the large-scale 
density fiuctuations. 

Fig. 9 shows the stellar mass function obtained by con- 
volving the model with a Gaussian whose mean and width 
varies with stellar mass according to the lines in Fig. 7. The 
fact that this matches very well the stellar mass function in- 
ferred from the rp-band mock catalogue demonstrates that 
this is the uncertainty that causes the systematic bias in the 
recovery of the stellar mass function. The J-band result is 
very similar. 

From the mock catalogue, it is apparent that a SWML 
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Figure 8. Systematic errors in the recovery of the global stellar 
mass function. The solid line traces the stellar mass function in 
the model. The dotted line shows a SWML estimate of the stel- 
laj: mass function from the mock catalogue using the true stellar 
masses, rather than those inferred from the galaxy light. Long- 
and short-dashed lines trace the normalised SWML stellar mass 
functions estimated from the mock, using the J- and rp-band- 
inferred stellar masses. The abundance drops for galaxies with 
log JO (Mstars/ h~-^ Mq) < 9 because of the blue luminosity limit 
in the model galaxy population, M(,j < —16. 

estimator of the stellar mass function is preferable to 1 /Knax, 
and that the large uncertainties in the stellar mass esti- 
mation lead to an overestimate of Qatara/i by ~ 40 per 
cent. Given that this overestimation results largely from the 
method chosen to assign stellar masses, it seems likely that 
such a systematic error would also be present when this tech- 
nique is applied to real data. While the size of this error is 
merely comparable with the systematic uncertainty arising 
from the choice of stellar IMF, it is larger than the statisti- 
cal uncertainties quoted by COl and Bell et al. (2003), and 
should be corrected for. Note that Bell et al. use a different 
procedure to calculate stellar masses, so the bias introduced 
by their method may not be quite the same as that found 
in this work. 

6.2 Comparison with other studies 

Applying the 1/Knax and SWML estimators to the full 
2dFGRS-2MASS sample within z = 0.12 yields the results 
shown in Fig. 10. Comparing with the results of COl, who 
used both Kennicutt (1983) and Salpeter (1955) IMFs, it 
is apparent that these new results are concordant at higher 
stellar masses. For the smaller systems, the J-band mass 
function, which should be directly comparable with that of 
COl, shows a deficit. The rp-band function does include 
most of the low stellar mass systems found by COl. Note 
that the stellar masses of the best-fitting Schechter function 
of COl have been multiplied by 1.1 to take into account. 
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Table 3. Best-fitting Schechter function parameters for the different galaxy SWML stellar mass functions. The SWML estimates are 
normalised using a x fit to the l/Vmax-estimatcd abundances over the range 10.4 < logio(Mstars//!.~^M0) < 11.4. A fit is then 
performed to determine the best-fitting Schechter functions. This latter fit is carried out over larger ranges of stellar masses. For the J- 
inferred case, galaxies within 10.0 < log^g (Avatars/ h~'^ Mq) < 11.4 are used, and this is extended to 9.5 < log]^()(Mstars/ Mq) < 11.4 
for the rp-inferred function. Also listed axe the mean Universal stellar mass densities derived from the SWML functions. Only statistical 
uncertainties are quoted. More important systematic uncertainties are discussed in the text. The final column contains the corrected 
estimates of the mean stellar density, talking into account the overestimation described in Section 5. Note that these stellar masses are 
all derived under the assumption that a Kennicutt IMF is universally applicable. 
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Figure 9. How stellar mass uncertainties affect the estimated 
stellar mass function. The solid line traces the stellar mass func- 
tion in the model. The dashed line shows the SWML estimate 
of the stellar mass function from the mock catalogue using the 

rp-band stellar masses. Points show the effect of convolving the 
model stellar mass function with a mass-dependent Gaussian, of 
mean and width given by the straight lines in Fig. 7. 

roughly, the difference between Kron and total magnitudes. 
The discrepancy with COl is mainly for the low-mass galax- 
ies, and thus does not greatly impact upon the inferred value 
of fistars/i. Assuming a Kennicutt IMF, COl found Slstars'i = 
(1.6 ± 0.24) X 10^3, whereas the corresponding numbers 
found here are (1.48 ±0.01) x 10"^ and (1.31 ±0.04) x 10"^ 
for the rp- and J-band respectively. As described above, it 
is entirely conceivable that those, already lower, estimates 
of fistarsfe are nevertheless still overestimated by ~ 40 per 
cent as a result of the errors inherent in assigning stellar 
masses to galaxies. Correcting these estimates under the as- 
sumption that they suffer the same fractional overestimates 
as the mock data, gives the corrected flstars/i values in the 
final column of Table 3. Also listed in Table 3 are the best- 
fitting Schechter function parameters obtained by a At to 
the normalised SWML stellar mass functions over the range 



Table 4. The recovered stellar mass functions from the real data, 
uncorrected for the measurement bias described in the text. The 
units of stellar mass arc h^^ Mq and the quoted abundances are 
dn/dlog]^Q(Mstars) in (h^^Mpc)^'^ . To correct for the bias intro- 
duced by the uncertainties in estimating stellar masses, applying 
a global shift of —0.15 in logjo(Mstars) is a good approximation. 
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9.5 < logio(Mstar3//i"^MQ) < 11.4 for the rp band stellar 
masses and 10.0 < log^Q(Mstars/ ^i"^ M©) < 11.4 for the J 
band stellar masses. The quoted statistical uncertainties are 
estimated from the scatter in the fits to the jack-knife sub- 
samples. In addition to the systematic uncertainties in the 
luminosity function estimation, a further source of system- 
atic uncertainty is the poorly known stellar IMF. To give 
a rough idea of how large an effect this can have, COl find 
that a Salpeter IMF yields a mean stellar mass density that 
is almost twice as large as that implied by the Kennicutt 
IMF. 

The difference between the stellar mass functions in- 
ferred using the two different mass estimators is consistent 
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Figure 10. The uncorrected stellar mass functions calculated 
using the real 2dFGRS-2MASS sample of galaxies. SWML re- 
sults from both J- (filled squares) and rp-band (black line with 
errors) stellar mass estimators are shown. Also shown are the 
best-fitting Schechter functions advocated by Panter, Ifeavens & 
Jimenez (Salpeter IMF; green); COl (Salpeter; solid blue); COl 
(Kennicutt; dotted blue); and Bell et al. - /f-band (solid red) and 
g-band (dotted red). 

with what was found by Bell et al. - namely that the sample 
selected by 2MASS flux is incomplete as a result of missing 
low surface brightness galaxies. This was shown by Bell et 
al. comparing a 2MASS if-selected sample with a SDSS g- 
band one. Their results exhibit a deficit of 2MASS-selected 
galaxies at log^Q[Mstars/(/i~^Af0)] < 10, very similar to the 
mass at which the J- and rp-band results found here begin 
to differ. 

Combining the estimate of the mean stellar density and 
the luminosity density from the previous section yields an 
estimate of the mean stellar mass-to-light ratio. Taking the 
average of the uncorrected r2stars(7'p) and f2stars(>/) values, 
because it is unclear which is better, gives a mean Ks band 
stellar mass-to-light ratio of 0.55Tq. If the corrected stellar 
mass densities are employed instead, then this reduces to 
0.39Tq. Even the uncorrected value is ~ 30 per cent lower 
than that found by COl (0.73 ± O.ISTq). 



7 STELLAR MASSES IN DIFFERENT SIZED 
GROUPS 

The 2PIGG catalogue allows a detailed study of how the 
galaxy stellar mass is distributed among groups of different 
size. This Section presents measurements of the stellar mass- 
to- light ratio and the stellar mass fraction, both as functions 
of group size. 

When calculating total group stellar masses, it is nec- 
essary to include that contributed by group members that 
do not make it into the flux-limited sample. This is done 
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Figure 11. Stellar mass-to-light ratios (Ks band) in the groups. 
The horizontal line shows the mean value in the model, and 
the solid curve shows the variation of the median with group 
size present in the model. Dashed and dotted lines trace the re- 
sults recovered from the mock catalogues using rp- and J-inferred 
stellar masses respectively. The corresponding median values in 
the 2PIGG sample are shown with filled squares and crosses. 
Taking the median group total mass-to-light ratio of Eke et al. 
(2004b) maps the following luminosities logjQ [L^j /(/i~^L0)] = 
(10, 11, 12) to masses logip[MQ/(/i-iMQ)] = (12.0, 13.6, 14.7). 

through the global bivariate SWML distribution of stellar 
mass and galaxy luminosity. Integrating this distribution 
over all stellar masses gives the fraction of stellar mass above 
a particular galaxy luminosity. The reciprocal of this frac- 
tion is the factor by which the total observed stellar mass 
must be increased in order to obtain the total group stellar 
mass. 

7.1 Group stellar mass-to-light ratios 

The stellar mass-to-light ratio in the Ks band is shown as a 
function of total group fo,j luminosity in Fig. 11. Only groups 
with 2 < Zmax = 0.05 + 0.02 [ logig (^-i),, / /i"^ L© ) - 10] are 
used, this being an appropriate compromise between having 
enough groups, and being able to see most of their stellar 
mass and luminosity. It is apparent that the values recov- 
ered from the mock generally overestimate the model mass- 
to-light ratio, particularly for the J-inferred stellar masses. 
This can be traced directly to the results shown in Figure 7. 
The slightly greater overestimation in stellar masses at the 
high mass end for the J-band estimates is sufficient to ac- 
count for this difference. A steeper trend is seen in the rp- 
inferred stellar mass-to-light ratio than in the J-band one. 
This is because for the less luminous groups, the rp mea- 
surements include groups with no galaxies that are detected 
in the near infra-red, and these groups preferentially have 
their stellar masses underestimated, thus dragging down the 
median rp-inferred stellar mass-to-light ratios. 
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Figure 12. The ratio of stellar to total dynamical mass inferred 
using the rp (upper panel) and J (lower panel) band estimators. 
A horizontal line shows the ratio Ostars /^m in the model, and the 
solid curve traces the variation of the median stellar mass fraction 
with group size in the model according to the semi-analytical 
prediction. A dotted line follows the results recovered from the 
mock catalogue, and the points show the median behaviour found 
in the 2PIGG sample. The dashed line traces the results found 
by LMS03. 



Very similar trends are apparent in the real 2PIGG 
data, with a small offset to lower mass-to-light ratios rel- 
ative to the mock. For the largest 2PIGGs, the stellar mass- 
to- light ratio is ~ O.GST©. However, if one assumes that 
this is overestimated by the same factor as in the mock cat- 
alogue, then one should expect that the underlying stellar 
mass-to-light ratio in the most luminous 2PIGGs should be 
~ O.45T0. 

Even with the ovcrcstimation of the stellar masses, the 
inferred stellar mass-to-light ratios are still lower than was 
assumed by LMS03 (~ O.8T0 for the richest clusters). This 
value was derived from dynamical data coupled to a maxi- 
mum stellar mass model, so it represents an upper limit on 
the stellar mass-to-light ratio. Alternatively, if the Kenni- 
cutt IMF were changed to one producing more stellar mass 
at a given luminosity (for instance the 'diet' Salpeter IMF 
of Bell & de Jong 2001), then the stellar mass-to-light ratio 
found here could be increased substantially to agree with 
that assumed by LMS03. 



7.2 Group stellar mass fractions 

Fig. 12 shows how the group stellar mass fractions (i.e. 
the mass currently in stars divided by the total group 
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Figure 13. Corrected stellar mass fractions in the real data (sym- 
bols), compared with the model (solid line). The dotted line is a 
least-squares fit to the stars), and the dashed line shows the re- 
sults of LMS03. 



mass) depend on group size for the rp- and J-inferred 

stellar masses. The horizontal line shows the mean stel- 
lar mass fraction in the model universe, and the other 
solid line traces the variation in stellar mass fraction with 
group size in the model, parametrised through the total 
group bj band luminosity. The median values recovered 
from the mock catalogues are shown with a dotted line. 
In constructing this curve, only groups with z < Zuiax = 
0.02 -I- O.O4[logio(i6j//i"^L0) - 10] are used. About 100 
groups contribute to the bin with logi,j[L6j/( h^^ L©)] = 11 
and ~ 25 groups to the bin at logj^Q[Lj,,/( h~'^ L©)] = 11.9. 
For groups with \ogiQ[Lb-^ /{h^^ Lq)] > 10.7 the variation 
of stellar mass fraction is recovered well, albeit with a ~ 25 
per cent offset resulting largely from the overestimation of 
stellar masses. This trend of increasing stellar mass fraction 
with decreasing group size is apparent in both the model 
and the mock catalogue. The 2PIGG data points show a 
generally similar behaviour to that in the mock catalogue. 
As for the case where stellar masses are inferred from the 
rp band, the J-band stellar mass estimator yields a trend of 
increasing stellar mass fraction with decreasing group size 
for the data. 

The dashed lines in Fig. 12 show the results of LMS03. 
They measure the mass fraction within a smaller radius than 

is used to define groups here. Thus, to calculate where their 
results fit on these plots, it has been necessary to assume 
a conversion from their M^oo to virial mass, Afvir, (this is 
taken simply as a factor 1.9 for a typical cluster dark matter 
density profile with a concentration, ca = 5), to use the 
typical mass-to-6j band light as a function of group size 
found by Eke et al. (2004b), and to account for their use of 
hYo- For the largest clusters, they find a stellar mass fraction 
of ~ 9 X 10"^/i"\ in contrast to the ~ 6.3 x 10"^/i"^ found 
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here. This is similar to the difference between the stellar 
mass-to-light ratios in the two studies. 

Interestingly, the variation of stellar mass fraction with 
cluster size found by LMS03 exceeds that found here, despite 
their clusters being in a region where no significant trend is 
apparent in the 2PIGG results. This is reminiscent of the sit- 
uation with the total mass-to-light ratio comparison. In fact, 
the 2PIGG results would show only weak evidence of any 
trend if it were not for the data at logj^Q[Lj,,/( L©)] — 
10.7. One could envisage that a strong trend of decreasing 
stellar mass fraction with increasing total mass could be 
artificially created if the lowest mass clusters preferentially 
had their masses underestimated. 

Under the assumption that the real mass fractions arc 
overestimated in the same way as those in the mock cat- 
alogue, a corrected stellar mass fraction can be derived. 
Fig. 13 shows a comparison between the model from which 
the mock was created and the corrected results for the real 
data. The different symbols show the two different stellar 
mass estimates and the mean (in log space) of these two. A 
dashed line traces the least-squares fit to the mean values. 
This line has the equation: 

logio(Mstars/i/-^) = -0.38 - 0.16 logio[i6j/( h-^ L0)].(7.1) 
For the biggest clusters, the typical corrected stellar mass 



fraction is only ■ 
the raw value. 



5 X 10 /i , about 25 per cent lower than 



8 WHERE ARE THE STARS? 

Having determined both the total stellar mass function of 
galaxies and the variation of stellar mass content with group 
size, the next natural step is to determine how all stars are 
partitioned among groups. It should be noted that this work 
ordy considers stars that arc associated with galaxies, ignor- 
ing any population of intergalactic stars. In large galgixy 
clusters, a few tens of per cent of the total cluster stellar 
mass may reside in intergalactic stars (e.g. Zibetti & White 
2004; Lin & Mohr 2004). However, only a small fraction of 
all stars reside in such large clusters, as will be shown. 

To determine the fraction of stars in groups of a given 
luminosity, a trivariate SWML method has been used. By 
choosing group luminosity to describe the group size, rather 
than mass, it is possible to use groups containing only a 
single galaxy. These contribute a significant fraction of the 
galaxies, and hence stellar mass, and are thus crucial in 
forming a full picture of where all of the stars reside. The 
trivariate SWML method extends the method used in the 
previous sections by adding an extra variable, namely the lu- 
minosity of the group in which any particular galaxy exists. 
Following the nomenclature of Loveday (2000), the space 
density of galaxies, <f>{Lg,L,M), with luminosity L, stellar 
mass M and residing in a group with total 6j luminosity 
Lg, can be used to write the probability of galaxy i taking 
particular values as 



Pi = 



<t>{LlL\M') 



g,min min 

The lower integration limit for the group total 6j luminosity 
takes account of the fraction of group luminosity that lies 



in galaxies beneath the flux limit of the 2dFGRS, assum- 
ing that group galaxies follow the Scliechter function deter- 
mined by Norberg et al. (2002) for the entire population 
(as was done by Eke et al. 2004b when defining group lu- 
minosity). The galaxy luminosity and stellar mass variables 
are either 6j-band luminosity and rp-band stellar mass or 
J-band luminosity and stellar mass. As shown by the inte- 
gration limits, all stellar masses are accessible at any given 
galaxy luminosity - this was also assumed for the bivariate 
case earlier. Performing the usual operations to maximise 
this probability while describing the galaxy space density 
in a stepwise manner leads to the following estimate of the 
density of galaxies residing in groups with 



ALg/2 < Lg < + ALg/2, 
having luminosity 

- AL/2 <L<L^ + AL/2, 



and stellar mass 

- AM/2 < M < M' 

0jk( = 



AM/2 
rijki 



EN 
i=l 



,AL„ALAM 



E-'ig Y^«t Y^«M 
P=l Z_-/q=l Z_-/r = l 



rHip 



,AL~ALAM 



(8.2) 

(8.3) 

(8.4) 
.(8.5) 



rijk! represents the total rmmbcr of galaxies in bin (j, k, I) and 
H ipqr is a trivariate version of the ramp function as described 
by Loveday (2000). NLg, Nl and Nm are the numbers of 
bins in group luminosity, galaxy luminosity and galaxy stel- 
lar mass respectively. 

The trivariate function can be projected to recover the 
galaxy luminosity function or the stellar mass function, for 
example. However, by projecting along the individual galaxy 
luminosity and stellar mass directions, the distribution of 
stars among groups of different luminosity can be retrieved. 

Fig. 14 shows both the differential and cumulative dis- 
tributions of stellar mass among different sized groups. 
While the distributions recovered from the mock catalogue 
are quite similar to that present in the model, there is a 
trend in the mock of placing a higher fraction of the stars in 
larger groups as the upper redshift of the sample is increased. 
This is particularly apparent for the rp-band results, because 
the mock J-band sample is also restricted by the relatively 
high 2MASS flux limit. For the z < 0.05 samples, a lack of 
large local clusters in the mock catalogue means that the 
results significantly underestimate the fraction of all stars 
in groups with \ogiQ{Li,.) > 11. As the upper redshift limit 
is increased, the volume encloses a fairer sample of the mock 
universe. Furthermore, the correction for missing group lu- 
minosity increases, which leads to an overestimate of the 
abundance of the more luminous groups. Consequently, the 
stellar mass is too frequently assigned to larger groups than 
should be the case. To correct for these systematic biases, 
one can choose an upper redshift that varies with group lu- 
minosity. The points in Fig. 14 show how the redshift limit 
can be chosen so that the model distribution is very ac- 
curately recovered. When the cut on group luminosity at a 
given redshift is introduced, it is necessary to alter the lower 
limit for the group luminosity integration in eq. (8.1). 

Assuming that similar systematic biases are appropriate 
in both the mock and real data, the points in the bottom row 
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Figure 14. The fraction of all stellar mass within groups as a function of group size, parametrized by group bj-band luminosity. 
Differential and the corresponding cumulative distributions are shown in the left and right columns respectively. The solid line in each 

panel shows the model distribution, while the other curves in the top two rows trace the estimates recovered as the maximum rcdsliift is 
changed from 0.05 (dotted) to 0.07 (short dashed) and 0.12 (long dashed). Results of using rp-infcrrcd stellar masses in the mock catalogue 
are shown in the top row, whereas the middle row contains the ./-inferred stellar masses. The symbols show the estimates with a variable 
upper redshift limit, dependent on the total group 6j band luminosity, as quantified in the legend. The bottom row shows the 2PIGG 
results, with the variable upper redshift limits for both the rp- and J-inferred stellar mass estimates, as well as tlic reference model. 
Taking the median group total mass-to-light ratio of Eke et al. (2004b) maps the following luminosities \ogiQ[Lf,-j / {h~ ^ Lq)] = (10, 11, 12) 
to masses logio[M0/(;i-iM"0)] = (12.0, 13.6, 14.7). 
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of Fig. 14 should represent unbiased estimates of the distri- 
bution of stellar mass among 2PIGG groups. The rp-band 
curves are very similar to those recovered from the mock 
catalogue, with the main differences being that the real Uni- 
verse contains a higher fraction of stars in very small groups 
and groups with \og^Q{L),-) ~ 10. The model stars are more 
frequently placed in groups with logjo ) ~ 9 — 9.5 than 
is the case in the real Universe. The difference for the small- 
est groups can be understood by recalling that the mock 
catalogue did not include galaxies with Mtj > —16. With 
these very low luminosity galaxies missing, fewer stars are 
placed into the smallest groups. In the J-band, the results 
arc broadly similar to the rp band, with a slight deficit of 
stars in small groups. This reflects the difference in faint-end 
slopes of the rp-inferred and J-inferred stellar mass func- 
tions shown in Fig. 10, at least some of which results from 
the fact that 2MASS misses some low surface brightness 
galaxies (Andreon 2002; Boll et al. 2003). As this additional 
systematic bias is known to affect the J-band results, the 
following fit is provided only for the rp-inferred stellar mass 
distribution: 

/(M,tar.)(> Lbj) = 1 - / dx, (8.6) 

where a = 0.96 and x = 10.39. This fit to the fraction of 
stellar mass iu haloes with a total blue luminosity exceeding 
is accurate to better than 0.02 for all group luminosities. 



up in stars. In conjunction with knowledge about what frac- 
tion of stellar mass is recycled, this gives a constraint on the 
integral over time of the Universal star formation rate. 

The stellar mass-to-light ratio and stellar mass fraction 
arc both studied as functions of 2PIGG size. Rich clusters 
are found to have stellar mass-to-ifs band light ratios of 
~ 0.6Tq , a value that is about 60 per cent larger than that 
found to be typical of groups with logjQ[L6j/(/i~^L©)] = 10. 
Removing the systematic errors inherent in this determina- 
tion reduces the cluster value to O.45T0. The opposite trend 
is found for the stellar mass fraction, with larger groups 
having smaller values. Taking into account the systematic 
error in the recovery gives a typical stellar mass fraction of 
~ 5 X W^'^h for the richest clusters. 

Finally, a trivariate stepwise maximum likelihood 
method is employed to paxtition, into groups of different 
size, the stellar mass that resides in galaxies in the local Uni- 
verse. It is found that only a couple of per cent of this stellar 
mass resides in galaxy clusters with log^Q[Z/(,j/( L©)] > 
12 (Iogio[M/(ft-^M0)] > 14.7). Half of this stellar mass 
is located in groups with logjQ[L6j/( Lq)] > 10.4 
(logjQ[M/( h^^ Mq)] > 12.5). Once again, using mock cata- 
logues can reduce the systematic biases associated with the 
entire measurement procedure. This adds weight to the as- 
sertion that the measured distribution is an accurate repre- 
sentation of the true underlying distribution of stellar mass 
among different sized haloes and, as such, represents a valu- 
able link in the chain connecting theory with observation. 



9 CONCLUSIONS 

The near-IR light and stellar content of the Universe as a 
whole, and its constituent groups, have been quantified us- 
ing the 2dFGRS, the 2PIGG catalogue and the 2MASS. This 
work extends that of previous studies with the largest sets 
of galaxies and groups yet used for these purposes. Further- 
more, the mock galaxy catalogues that have been employed, 
allow a careful quantification of the systematic errors as- 
sociated with estimating stellar masses. This represents an 
important advance, because the systematic bias introduced 
by errors in the estimated stellar mgisses, which has been 
quantified via the mock catalogues, is sufficiently large that 
the mean stellar mass density had previously been signifi- 
cantly overestimated. The mock catalogues permit this bias 
to be quantified and corrected for. 

The mean luminosity density in the Universe is found 
to be pj = (3.57 ± 0.11) x lO^ftLgMpc"^ and pxs = 
(7.04 ±0.23) X lO^/iLoMpc"^ (statistical uncertainty only) 
in the two near-IR bands considered here. Systematic uncer- 
tainties are likely to be of order 20 per cent and dominated 
by the completeness of the catalogues used and the con- 
version of Kron to total magnitudes. Taking into account 
the ~ 40 per cent overestimation of the stellar mass den- 
sity that, according to the mock catalogues, arises from the 
uncertainty in inferring stellar masses from galaxy fluxes 
and colours, the mean stellar mass density amounts to 
fistars/i = (0.99 ± 0.03) X 10~^. This value assumes that a 
Kennicutt IMF is applicable to all sites of star formation. If 
a Salpeter IMF had been chosen then the estimated fistars/i 
would be a factor of ~ 2.1 larger. The value of fistars/i pro- 
vides a measure of the density of material currently locked 
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